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ABSTRACT 

The first galaxies that formed at a redshift z ~ 20-30 emitted continuum photons with energies between the 
Lya and Lyman limit wavelengths of hydrogen, to which the neutral universe was transparent except at the 
Lyman-series resonances. As these photons redshifted or scattered into the Lya resonance they coupled the 
spin temperature of the 21cm transition of hydrogen to the gas temperature, allowing it to deviate from the 
microwave background temperature. We show that the fluctuations in the radiation emitted by the first galaxies 
produced strong fluctuations in the 21cm flux before the Lya coupling became saturated. The fluctuations were 
caused by biased inhomogeneities in the density of galaxies, along with Poisson fluctuations in the number of 
galaxies. Observing the power-spectra of these two sources would probe the number density of the earliest 
galaxies and the typical mass of their host dark matter halos. The enhanced amplitude of the 21cm fluctuations 
from the era of Lya coupling improves considerably the practical prospects for their detection. 

Subject headings: galaxies: high-redshift, cosmology: theory, galaxies: formation 



1. INTRODUCTION 

The reionization history of cosmic hydrogen, left over 
from the big bang, provides crucial fossil evidence for when 
the first stars and black holes formed in the infant universe 
(Barkana & Loebl 120011) . The hyperfine spin-flip transition 
of neutral hydrogen (H I) at a wavelength of 21 cm is po- 
tentially the most promising probe of the cosmic gas before 
reionization ended. Observations of this line at a wavelength 
of 21 x (1+z) cm can be used to slice the universe as a func- 
tion of redshift z and obtai n a three -dimensional map of the 
diffuse H I distribution in it (Hogan & Rees 1979). The 21cm 
signal vanished at redshifts z > 200, when the small resid- 
ual fraction of free electrons after cosmological recombina- 
tion kept the gas kinetic temperature, Tj, close to the temper- 
ature of the cosmic microwave background (CMB), T 1 . Sub- 
sequently, between 200 > z > 30 the gas cooled adiabatically, 
faster than the CMB, and atomic collisions kept the spin tem- 
perature T s of the hyperfine level population below r 7 , so that 
the gas appeared in absorption (IScott & Reeslll990l) . Primor- 
dial density inhomogeneities imprinted a three-dimensional 
power-spectrum of 21cm flux fluctuations on scales down to 
5 10 comovingkpc, making it the richest data set on the sky 
(Loeb & Zaldarriaga 2004). As the Hubble expansion contin- 
ued to rarefy the gas, radiative coupling of T s to T 7 dominated 
over collisional coupling of T s to 7j and the 2 1 cm signal began 
to diminish. 

However, as soon as the first galaxies appeared, the 
UV photons they produced between the Lya and Ly- 
man limit wavelengths propagated freely through the uni- 
verse, redshifted or scattered into the Lya resonance, and 
couple d T s and 71 once a gain through the Wouthuysen- 
Field (IWouthuvsenl 119521 iFieldl 119581) effect by which 
the two hyperfine states are mixed t hrough the absorp- 
tion and re-emission of a Lya photon (iMadau et alJ 119971 
Ciardi & Madau 2003). Emission of UV photons above the 
Lyman limit by the same galaxies initiated the process of 
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reionization by creating ionized bubbles in the neutral gas 
around these galaxies. While UV photons were blocked at 
the neutral boundary of the bubbles, X-ray photons propa- 
gated farther into the bulk of the intergalactic gas and heated 
Tk above T 7 throughout the universe. Once T s grew larger than 
T 7 , the gas appeared in 21cm emission with a brightness tem- 
perature that was asymptotically independent of T s and of 7j. 
The ionized bubbles imprinte d a knee in the power-spectrum 
of 21cm fluctuations (Zaldarriaga et al. 2004), which traced 
the topology of H I as the ionized bubbles percolated and com- 
pleted the process of reionization dFurlanetto et al.l20 04a). 

Since the first galaxies represented rare peaks in the cos- 
mic density field, their spatial distribut ion possessed unusu - 
ally large fluctuations on large scales (Barkana & Loeb 2004). 
In this paper we show that the resulting large-scale bias of the 
first galaxies produced substantial 21cm fluctuations during 
the intermediate epoch when the first Lya photons coupled T s 
to Tk- Since a relatively small number of galaxies contributed 
to the flux seen at any given point, Poisson fluctuations were 
significant as well, and we show that they produced correlated 
21cm fluctuations. Owing to the rapid evolution in the den- 
sity of galaxies with cosmic time, there is a sharp decline in 
the abundance of sources with distance (corresponding to re- 
tarded cosmic time for the source emission due to the light 
travel time) around each point in the intergalactic medium. 
Ultimately, there is an effective horizon around each point out 
to which sources may contribute to 21cm fluctuations. In ad- 
dition, the atomic physics of resonant scattering by hydrogen 
fixes different horizon distances for photons emitted in vari- 
ous wavelength intervals between Lya and the Lyman limit; 
these horizons are set by the requirement that a given photon 
must redshift or scatter into the Lya resonance at the point 
of interest. We incorporate all of these ingredients into our 
calculation. 

The outline of the paper is as follows. In § 2 we summa- 
rize the basic equations that describe the evolution of the spin 
temperature. In § 3 we review the angular anisotropy of the 
21cm power spectrum due to peculiar veloci ties, as derived in 
our previous paper (Barkana & Loeb 2005). We then calcu- 
late the flux of Lya photons produced by the earliest galaxies 
(§4) and its fluctuations (§5). Finally, we demonstrate how 
the resulting power spectrum of fluctuations in the 21cm flux 
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(§ 6) can be used to study the properties of the earliest galaxies 
(§ 7), and briefly summarize our main results (§ 8). Through- 
out the pap er we assume the co ncordance set of cosmological 
parameters (Spergel et al. 2003). 

2. SPIN TEMPERATURE HISTORY 

The H I spin temperature, T s , is defined through the ratio 
between the number densities of hydrogen atoms in the ex- 
cited and ground state levels, n\/no = (gi/gu)^p{~T*/T s }, 
where (g\/go) = 3 is the ratio of the spin degeneracy fac- 
tors of the levels, and T+ = 0.0682K. When the spin tem- 
perature is smaller than the CMB temperature, neutral hy- 
drogen atoms absorb CMB photons. The resonant 21cm ab- 
sorption reduces the brightnes s temper ature of the CMB by 
dScott & Reesll 990UMadau et alF997) 



1+z 



(1) 



where the optical depth for resonant A = 21cm absorption is, 

3cA 2 Mio«h 
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where tin is the number density of hydrogen, A\q = 2.85 x 
10~ 15 s _I is the spontaneous emission coefficient, and dv r /dr 
is the gradient of the radial velocity along the line of sight 
with v r being the physical radial velocity and r the comov- 
ing distance; on average dv r /dr = H(z)/(l +z) where H is the 
Hubble parameter. The velocity gradient term arises since the 
21cm scattering cross-section has a fixed thermal width which 
translates through the redshift factor (1 + v r /c) to a fixed inter- 
val in velocity (Sobolev 1960). In this paper we consider the 
epoch long before complete reionization, and we assume that 
the hydrogen gas is almost entirely neutral. 

The mean brightness temperature offset on the sky at red- 
shift z is 



7i = 28mK 
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(3) 

where we have substituted the concordance values for the cos- 
mological parameters 51;,, h, and fi m . The spin temperature 
itself is coupled to 71 through the spin-flip transition, which 
can be excited by atomic collisions or by the absorption of 
Lya photons. Assuming that the photons near Lya are redis- 
tributed in frequency according to the thermal velocities of the 
atoms (i.e., the spectrum is thermalized according to 71), the 
combination that appears in T/, becomes (Field 1958) 



T s -T„ 



x tot 
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where the sum of the radiative and collisional 

threshold parameters. These parameters are 
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where P a is the Lya scattering rate which is proportional to 
the Lya intensity, and Avi_n(71) is tabulated as a function of Tj. 
(Allison & Dal garnolH969t IZvgelmanl 12005). Note that we 
have adopted a modified notation compared to the y a and 



y c previously used (lFieldll958t iMadau etal1ll997l) . defining 
x a = y Q (r Q /r 7 ) and x c = y c (7^/r 7 ). The coupling of the spin 
temperature to the gas temperature becomes substantial when 
Xt n t > 1; in particul ar, x a = 1 defines the thermalization rate 
(Mada u et al.lll997l) of P a , which is 7.69 x 10" 12 (1 +z)/10 
photons per second. 

3. THE SEPARATION OF POWERS 

Although the mean 21cm emission or absorption is difficult 
to measure due to bright foregrounds, the unique character 
of the fluctuati ons in 71 allows for a much easier extraction 
of the signal (iGne din & Shaverl20 04l:IZa ldarriaga et all2004t 
iMorales & Hewittll2004t iMorale si 12004 ISantos et alJl2004l) 
During the era of initial Lya coupling, fluctuations in Tb are 
sourced by fluctuations in the density {5), temperature (St,,) 
and radial velocity gradient of the gas as well as the Lya flux 
through x a (5 Xa ). The fluctuations in the neutral fraction of 
the gas become important only later, when the volume filling 
fraction of ionized bubbles is substantial. 

The above sources of fluctuations are all isotropic except 
for the velocity gradient term. Since this term involves ra- 
dial project ions, it inserts an anisotropy into the power spec - 
trum ( Kais er! 1987tlTozzi et al.l2 000: Bharadwai & Ali 2004). 
The to tal power spectrum can be written as (Barkana & Loeb 
2005) 



P Th (k) = n 4 P s (k) + 2fi 2 



pP s (k) + ^P S - a (k) 

Xtot 



x tot 
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P a (k)+^Ps- a (k) 
Xtot 
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where i tot = (1 +A:tot)^tot, Ps and P a are the power spectra of 
the fluctuations in density and in x a , respectively, and the 
power spectrum P$- a is the Fourier transform of the cross- 
correlation function, 



Here (3 is given by ( Barka na" & Loebl2005h 
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where the adiabatic index is 7,, = 1 + (StJS). Also, /1 is 
the cosine of the angle between the wavenumber k of the 
Fourier mode and the line-of-sight direction. Considering 
equation Q as a polynomial in [i, i.e., 



P Tb (K) = y A P^ + H 2 Pf+Pn° , 



(10) 



we note that the power at just three values of /1 is required in 
order to separate out observationally the coefficients of 1, fi 2 , 
and /i 4 for each k. In particular, we can construct a combina- 
tion that probes whether some sources of flux fluct uations are 
uncorrected with S, through ( Bark ana & Loebl2005l) 
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Note that since 21cm measurements on large angular scale s 
cannot yield significant information (Barkana & Loeb 2005), 
we focus in this paper on three-dimensional power spectra, 
which are directly measurable from 21cm observations in 
small angular fields. 

We note that several additional redshift effects do not sig- 
nificantly modify the 21cm brightness fluctuations. First, 
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consider an observed correlation function (or the equivalent 
power spectrum) that is calculated by averaging the prod- 
uct of the 21cm temperatures at pairs of points within a 
small three-dimensional volume. We assume that this vol- 
ume is centered on redshift z and has a redshift width Az 
that corresponds to the observed frequency range used to cal- 
culate the correlation function at z. The change in the cos- 
mic mean temperature Tb(z') with redshift z' within this slice 
would produce variations in temperature even in an homo- 
geneous expanding universe. However, these variations are 
coherent on the sky and have a very simple form. Denot- 
ing by / the line-of-sight component of r, where / = marks 
the center of the redshift slice and AZ its width, the cor- 
relation function ^(r) = (5^ (•" 1)^(1*1 +r)) is modified to 
((5 Tb (r\) + lidiTb)[5T b (ri+r) + (li + l)diTb]), where we denote 
the derivative [dT/,(z)/dl] at z by d{T[ 7 . Noting that (5 Th {v\)) 
and are both zero within the slice, we find that £r t (r) re- 
ceives an extra term, {d[Tb) 2 (l 2 ). Since (I 2 ) = (A/) 2 /12, the 
finite redshift width Az adds to the power spectrum Pt,,(^) a 
term 



by 



^:A ; (k)=-vr 3 



Az 



dT b (z) 



dz 
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Thus, this effect does not modify the power spectrum at k =/ 0. 

In addition, peculiar velocities change the apparent distance 
along the line of sight. Over a distance I, the unperturbed ve- 
locity v = aHl receives a fractional perturbation S v (l) which 
is of order the density perturbation 6 averaged on the scale 
I (e.g., S v = -S/3 for a spherical tophat perturbation out to 
radius /). The mean temperature T b is modified due to this 
displacement by ~ l5 v {l)diTb, which is smaller than 8 by 
~ alH jc which at z = 20 is Z/(2Gpc) -C 1. Also, to first or- 
der we can neglect the effect of the peculiar velocities on the 
perturbation 8j b itself. Finally, gravitational potential fluctua- 
tions are ~ 10~ 5 on the scale of the horizon and are suppressed 
relative to the density fluctuations by a factor of ~ (alH/c) 2 
on small scales. Note that 10~ 5 fluctuations are important in 
CMB observations, since at z ~ 1100 density fluctuations on 
small scales had not yet grown due to the radiation pressure 
of the photon-baryon fluid; in the redshift range that we con- 
sider for 21cm observations, however, density fluctuations are 
already at a level of ~ 1% on 100 Mpc scales (at z = 20), and 
it would be ve ry difficult to measure the ~ 10~ 4 fluctuations 
on Gpc scales (Barkana & Loeb 2005). 

4. THE Lya FLUX OF GALAXIES 

The intensity of Lya photons observed at a given redshift 
is due to photons that were originally emitted between the 
rest-frame wavelengths of Lya and the Lyman limit. Photons 
that were emitted below Ly/3 by some source simply redshift 
until they are seen by an atom at z at the Lya wavelength. 
Such photons can only be seen out to a distance corresponding 
to the redshift z ma x(2), where 1 +z ma x(2) = (32/27) x (1 + z), 
where 32/27 is the ratio of Lya to Ly/3 wavelengths. Pho- 
tons above the Ly/3 energy redshift until they reach the near- 
est atomic level n. The neutral intergalactic medium (IGM) is 
opaque even to the higher levels and so the photon is absorbed 
and re-emitted. During these scatterings, the photon is almost 
immediately downgraded to a Lya photon and then keeps on 
scattering, and so most of the atoms that interact with it en- 
counter it as a Lya photon. To be seen at the Lya resonance 
at z, the photon must have been emitted below a redshift given 



1 +Zmax(«) = (1 +Z) 
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Thus, we simply add up the number flux of photons emitted 
between consecutive atomic levels, and integrate over sources 
out to the distance corresponding to the appropriate redshift 
interval. 

Note that while sources are expected to also produce a 
strong Lya line, due to internal reprocessing of ionizing pho- 
tons into Lya photons, these photons are not expected to have 
a significant effect. The photons stream into the H II region 
of their source, and by the time they reach the edge, many 
of them have redshifted above the Lya wavelength by more 
than the line width; the remaining ones are trapped and dif- 
fuse only across a distance of ~ tens of parsecs before they 
too redshift out of the line. Thus, these photons only affect 
a very small fraction of the IGM, much less than the ionized 
fraction, which is itself small in the redshift range of inter- 
est here. Even though the Lya flux within this volume can 
be high, the brightness temperature saturates when x tot — > 00. 
Observations on scales larger than the size of an H II region 
(which is < 1 comoving Mpc) smooth the signal over many 
such regions. The overall correction is thus small both on the 
mean brightness temperat ure and on th e fluctuations in it. 

We therefore get (Mad auet all 19971) 



167r 2 7> 2 / a 
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where J a is the proper Lya photon intensity (defined as the 
spherical average of the number of photons hitting a gas ele- 
ment per unit a rea per unit time per unit frequ ency per stera- 
dian), and S a (IChen & M iralda-Escude 2004) is a correction 
factor between 1 and 2 that depends on the gas temperature 
and accounts for the redistribution of photon frequencies near 
the Lya line center due to scattering off the thermal distribu- 
tion of atoms. The intensity is itself given by a sum over n 
(see also §|6j, 
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where absorption at level n at redshift z corresponds to an 
emitted frequency at z' of 



^ = ^ll(1 



-n ) 

(1+z) 



(16) 



in terms of the Lyman limit frequency ^ll- The comoving 
photon emissivity (defined as the number of photons emitted 
per unit comoving volume, per proper time and frequency, at 
rest-frame frequency v at redshift z) is in turn 



e(u,z) = n b f* ^gai(z) £b(v) , 



(17) 



where if b is the cosmic mean baryon number density at z = 0, 
/* is the efficiency with which gas is turned to stars in galactic 
halos, Fgai(z) is the gas fraction in galaxies at z, and tbW) is 
the spectral distribution function of the sources (defined as the 
number of photons per unit frequency emitted at v per baryon 
in stars). 

To determine the gas fraction in galaxies, we must first 
find the halo mass functi on. Although the basic model of 
iPress & SchechteJ i ll 9741) describes the results of numeri- 
cal simulations qualitatively, we use the halo mass function 
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of LSheth & Tormer] (fl999h which fits numerical simulations 
more accurately Uenkins etafl feOOl). In order to calculate 
fluctuations among large-scale regions of various densities in 
the global gas fraction that resides in galaxies within them, 
we must also find the halo mass function as a function of the 
overall mean density in each region. We adjust the mean halo 
di stribution as a function o f density based on the prescription 
of iBarkana & Loebl J2004I) . where we showed that this pre- 
scription fits a broad range of simulations including ones run 
at z ~ 20. As gas falls into a dark matter halo, it can frag- 
ment into stars only if its virial temperature is above 10 4 K 
for cooling mediated by atomic transitions, or ~ 500 K for 
molecular H2 cooling. Given a minimum cooling mass, we 
consider two possibilities for the star formation efficiency: ei- 
ther a fixed efficiency in all halos, as in equation Mil , or an 
efficiency that drops as the square of the halo circular veloc- 
ity V r at V r 5, 180 km s" 1 , as su ggested by local observations 
( Deke l & Wool2002HKauffman et al 120031) : in the latter case, 
we integrate over the halo mass function in order to determine 
the net effective value of/*. 

The spectral distribution function €b(v) depends on the as- 
sumed sources of radiation. In our discussion, we assume 
that stars dominate over mini-quasars, as expected for t he low- 
mass g alaxies that are typical at high redshift ( Wvithe & Loeb 
2004b). The stellar spectrum depends on the initial mass func- 
tion (IMF) and the metallicity of the stars. In what follows, we 
co nsider two ex amples. The first is the locally-measured IMF 
of lScaiol ( 1998) with a metallicity of 1/20 of the solar value, 
which we refer to as Pop II stars. The second case, labeled as 
an extreme Pop III IMF, consists entirely of zero metallicity 
M > 100 Mm stars, as expected fo r the earliest generation of 
galaxies (Bromm & Larson 2004). In each case, we approx- 
imate the emissivity as a separate power law tbW) oc v as ~ l 
between every pair o f consecutive levels of atomic hydrogen. 
For example, using iLeitherer et alJ d 1 999 ') we find that the 
Pop II stars emit 9,690 photons per baryon between Lya and 
the Lyman limit, out of which 6,520 photons are between 
Lya and Ly/3 with an as = 0.14 in this lowest frequency in- 
terval. The cor responding numbe rs for Pop III stars [using 
iBromm. Kudritzki. & Loebl d200ll) l are 4,800 , 2,670 , and 
aj = 1.29. 

Figure [J illustrates the cosmic evolution of the various 
key temperatures. The residual electron fraction keeps the 
gas temperature 7* locked to the CMB temperature [T 7 = 
2.725 (1 +z) K] until l+z~ U2(VL h h 2 /0.024) 2 / 5 (iPeeblesI 
1993), but subsequently the gas cools adiabatically. The mean 
temperature offset 7j is ~ -30 mK at z = 100 and it rises to- 
wards zero at low redshift until the Lya flux from galaxies 
drops it again, to -100 mK and below. Note that this Figure 
and our results in subsequent sections assume an adiabatically 
cooling cosmic gas, except where we explicitly assume pre- 
heating by X-rays. 

5. FLUCTUATIONS IN Lya FLUX 

There are two separate sources of fluctuations in the 
Lya flux. The first is density inhomogeneities. Since grav- 
itational instability proceeds faster in overdense regions, the 
biased distribution of rare galactic halos fluc tuates much 
more than the global dark matter density (iKaiserl 119841 
Bark ana & Loe b 2004). When the number of sources seen 
by each atom is relatively small, Poisson fluctuations provide 
a second source of fluctuations, statistically independent of 
the first source, to linear order. Unlike typical Poisson noise, 
these fluctuations are correlated between gas elements at dif- 
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FIG. 1. — Redshift evolution of various mean temperatures. The mean 
21cm spin temperature is shown in the upper panel for adiabatic IGM gas 
with Lya coupling (solid curve) and without it (dotted curve). Also shown 
for comparison are the gas temperature (long-dashed curve) and the CMB 
temperature (short-dashed curve). The mean 21cm brightness temperature 
offset relative to the CMB is shown in the lower panel, with Lya coupling 
(solid curve) and without it (dotted curve). The Lya flux is calculated as- 
suming that galaxies with Pop III stars form in dark matter halos where the 
gas cools efficiently via atomic cooling. The star formation efficiency is nor- 
malized so that the Lya coupling transition (i.e., x tot = 1) occurs at redshift 
20. 



ferent places, since two nearby elements see many of the same 
sources, where each source is at about the same distance from 
the two gas elements. 

Since each hydrogen atom receives some Lya flux from 
sources as far away as ~ 250 comoving Mpc, the flux is 
naively expected to be highly uniform based on the homo- 
geneity of the universe on these large scales. However, the 
fluctuations in flux are actually relatively large because a sig- 
nificant portion of the flux comes from nearby sources. One 
reason is the 1 jr 2 decline of flux with distance. In an homo- 
geneous, non-expanding universe, the flux per log(r) would 
scale linearly with r, since the l/r 2 decline of the flux is over- 
come by the r 3 volume factor; while large scales would dom- 
inate the total flux, the relative contribution of small scales 
would already be far more significant than if we were sim- 
ply counting galaxy numbers, where the overall scaling would 
be cubic instead of linear in r. The l/r 2 scaling of flux also 
strongly magnifies Poisson noise since, e.g., the 21cm emis- 
sion along a given line of sight can fluctuate strongly due to a 
single Lya source that happens to lie very close to that line of 
sight. Although the mean flux at a point is well defined, if we 
calculate the mean squared flux at a point, then a shell at a dis- 
tance r contains a number of sources oc r 2 and contributes to 
the squared flux a value oc 1 /r 4 ; thus, the mean squared flux 
diverges at small radii, and since the correlation function of 
flux at r = is proportional to its mean square, it too diverges 
at small r. 

A second reason for the dominance of small scales in the 
Lya flux is that higher Lyman series photons, which are de- 
graded to Lya photons through scattering as discussed in the 
previous section, can only be seen from a small redshift inter- 
val that corresponds to the wavelength interval between two 
consecutive atomic levels of hydrogen. A third reason is that 
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distant sources are observed after a time delay determined by 
the speed of light, and their relative contribution is reduced 
since galaxy formation was not as advanced at earlier times. 
To facilitate the weighting of different scales, we summarize 
equations H4\ . (I15> . and (I17> in the form: 



dr dz' 



(18) 



The top panel of Figure [2] shows an estimate of the rela- 
tive contribution to the fluctuations in Lya flux from sources 
at various distances from a gas element at z = 20. We take 
the contribution to the flux per log(r), i.e., dx a /d\ogr (which 
is cx r at small distances), and multiply it by the mean den- 
sity fluctuation in spheres of radius r, in order to estimate the 
contribution per log(r) to the density-induced fluctuations in 
flux. We assume a scale-invariant spectrum of primordial den- 
sity fluctuations, as predicted by inflation and indicated by 
the latest observations (Selia ket alJl2005l) . Since this spec- 
trum yields fluctuations that diverge only logarithmically at 
small scales, the declining contribution of small scales to the 
flux dominates as r — ► 0. On the other hand, density fluctu- 
ations decrease as r~ 2 on scales r > 100 Mpc, which over- 
comes the cx r contribution to the flux and reduces the role 
of large-scale modes in the flux fluctuations. Although these 
effects make the ~ 10-100 Mpc scales dominate the density- 
induced fluctuations in Lya flux, the drop in the density of 
time-retarded sources (as evident from the comparison of the 
top two curves in the top panel), plus the large contribution of 
nearby sources through photons above Ly/3 , make the overall 
effect smoothed out over scales of ~ 0. 1-100 Mpc. The figure 
shows discontinuous jumps at distances that correspond to the 
various redshifts z max («) that mark the Lya horizons for pho- 
tons emitted with energies between the n — ► 1 and n + 1 — ► 1 
transitions of hydrogen. 

The bottom panel of Figure |2 shows the relative contribu- 
tion to Poisson-induced fluctuations in the Lya flux. We con- 
struct this estimate by multiplying dx a /d\ogr by the Poisson 
1 / \/N fluctuation, where for small distances N oc r 3 . In a 
non-expanding universe, this would already diverge as 1 / 1 %fr 
at small scales. The previously-mentioned effects only in- 
crease the dominance of small scales in the Poisson fluctu- 
ations. Thus, measuring the two different sources of flux fluc- 
tuations has a complementary effect and allows a broad range 
of scales to be probed. Furthermore, the Poisson fluctuations 
are independent (to linear order) of the underlying primordial 
density fluctuations, and more directly measure the number 
density of sources. In Figure|2]and the following figures, we 
include in the case of "full Lyman-band emission" all the pho- 
tons emitted from Lya up to a maximum range given by the 
interval between levels n max and n max + 1 of hydrogen (assum- 
ing that none of the photons are obscured or absorbed by dust 
within the host galaxy). We set this maximum level by cutting 
off photons that could only be seen from sources so nearby 
that the gas element would have to be within the H II region of 
the individual source, and thus this gas would be ionized and 
would not contribute to the 21cm optical depth. Calculating a 
typical H II region size based on an escape fraction of ionizing 
photons of / esc = 0.3, this yields a maximum level n max = 23 in 
the Pop III case at z = 20, and corresponds to including 99.3% 
of all photons emitted between Lya and the Lyman limit. We 
note that since ionizing sources are clustered JWyithe & Loebl 
2004b; Barkana & Loebl 12004 iFurlanetto et al.ll2004bl) . the 
appropriate cutoff scale might be somewhat larger. 
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FIG. 2. — Relative contribution per log(r) of shells at various radii r from 
a gas element at redshift z = 20, to the fluctuations in the Lyo flux. We con- 
sider the fluctuations sourced by density inhomogeneities (top panel) and by 
Poisson fluctuations (bottom panel). In each set of solid curves, we show 
from top to bottom (at large r) the asymptotic case of a non-expanding uni- 
verse followed by galaxies, where we include stellar radiation emitted up to 
Ly/3, Ly5, or the case of full Lyman-band emission (see text). The minimum 
galaxy mass is set by the atomic cooling threshold in dark matter halos. In- 
side galaxies, we assume a fixed star formation efficiency and a Pop III stellar 
population. Also shown for comparison is the case of full Lyman-band emis- 
sion with Pop II stars (dotted line). All curves are arbitrarily normalized to 
unity at r = 0. 1 Mpc. 



6. POWER SPECTRUM OF FLUCTUATIONS IN THE 21CM FLUX 

We first consider the effect of density perturbations on the 
flux as given by equation Jl 8i . When considering pertur- 
bations, we let r be the Lagrangian position, which corre- 
sponds to the comoving distance at some initial high red- 
shift. We now consider a field of perturbations <5o(r), lin- 
early extrapolated to the present from the initial time with 
the growth factor D(z'), so that the perturbations at a redshift 
z' are So(r)D{z'). Consider the contribution of a Lagrangian 
volume dV = dAdr at position r to the flux observed at the 
origin of the coordinates at redshift z- The emission rate 
e{v' ,z') is different than average, since the density perturba- 
tion So(r)D(z') changes the abundance of galaxies. On large 
scales, the change in the abundance of h alos is linear in the 
density perturbation, fo r every halo mass ( Mo & Whi te! 199H 
Barkana & Loeb 2004), yielding a weighted bias b(z') in the 
overall emission rate. Thus, the emission rate is increased by a 
factor [l+b(z')S (r)D(z')]. Note that the required emitted fre- 
quency is determined by the redshift just as in the unperturbed 
case [i.e., equation d!6b 1. and also to first order we can use 
the unperturbed z'(r) relation in the perturbation term because 
corrections to this are small. The observed proper intensity 
at z due to a proper emitting area dA p at the source at red- 
shift z' is unchanged by perturbations, due to the conservation 
of surface brightness known from gravitational lensing the- 
ory. However, for a fixed Lagrangian area dA, the proper area 
dA p (and thus also the observed intensity) is increased by the 
factor [1 + |(5()(r)D(z')] to linear order. Finally, we must relate 
the remaining factor in the volume element dV', i.e., the radial 
width dr, to the redshift interval dz'- The appropriate relation 
is dz' = [dv r /dr]dr/c, where v r is the radial component of the 
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physical velocity. The unperturbed form of this Doppler re- 
lation is dzf = H{z')dr/[(\ + z')c], which accounts for the de- 
nominator in equation J15I . Even in the perturbed case, the 
radial integration limit is fixed in redshift as in equation (I15> . 
Note that we have assumed that the retarded time correspond- 
ing to a given source redshift z' > z is unperturbed; the pertur- 
bation in peculiar velocity causes a Az ~ 8ar/(c/H), which 
is negligible compared to the biased gal axy fluctuation which 
corresponds to a Az ~ £(1 +z) (Barkana & Loeb 2004). 

Integrating over the volume dV', we add the various terms 
mentioned before and obtain the total linear perturbation in 
x a , which is a spatial average of the density perturbation field, 
plus a spatial average of the velocity gradient term. The frac- 
tional perturbation in dv r /dr has a simple Fourier transform, 
~5 drVr = -pi 2 5 JKaiserlll987t iBharadwai & Ali 2004-j), which is 
used to sepa rate the angular structu re of the power-spectrum 
[see §|3]and Barkana & Loeb (2005)]. Here, /i is the cosine 
of the angle between the wavenumber k of the Fourier mode 
and the r (i.e., line of sight) direction. After proper averaging 
we find the Fourier perturbation in x a in terms of the Fourier 
transform of the density perturbation at z. 

5 a (k) = W(Jt)S(k), (19) 

where 

W(k)=- f^ 2 \ z '^^Ul +b { z ^ m kr)-\i 2 {kr)\ , 
x a J, dz' D(z) y 3 J 

(20) 

expressed in terms of spherical Bessel functions. The first 
term is reminiscent of the effective bias in galaxy number den- 
sity, where the addition of 1 to the Lagrangian bias accounts 
for the higher physical density per unit volume. The second 
term is a geometrical correction for the case of flux, which 
arises from the velocity gradient term. The contribution to the 
power spectrum P a (k) is then W 2 (k)Ps(k), while the contribu- 
tion to P$- a {k) is W(k)Ps(k). We note that P a corresponds to 
a real-space correlation function of the flux at two points A 
and B separated by a comoving distance /. In this correlation, 
if point A is closer to sources at some location r then point B 
will see the progenitors of the galaxies seen by point A. The 
relation between the flux of the progenitors and the flux of the 
later parent is implicitly accounted for in the above weighting 
by the redshift-dependent bias factor b(z'). 

If the number of galaxies N that contribute to the flux is 
relatively small, then in addition to fluctuations in the mean 
number density and mass distribution of halos we expect sig- 
nificant Poisson fluctuations. We first derive a general frame- 
work for analyzing correlations induced by Poisson fluctua- 
tions. We note that in t he context of the ionizing intensity due 
to quasars, Zuo ( 1992) derived the joint probability distribu- 
tion of the intensities at two different points. Focusing on the 
two-point correlation function only, we present here a deriva- 
tion that is simpler and more general, including cosmological 
redshift and time-retardation effects. 

We consider a large number of small volumes dV , such that 
the mean number of halos (dn) <C 1 in each, so that dn 2 = dn 
since at most a single halo may be found within a single vol- 
ume element. Although we describe a volume distribution 
for convenience, this setup also applies to a generalized no- 
tion of volume, where the population of halos is divided up 
into many small units, and the effective position within this 
"volume" is specified by as many parameters as needed, in 
addition to spatial position (e.g., halo mass, emission time, 



merger history, etc.). We denote by dx^ the flux that ar- 
rives at a point A due to dn halos in some volume element; 
this flux is proportional to dn, with a proportionality factor 
that may depend on the distance to A and on the various pa- 
rameters that specify the volume element, and can incorpo- 
rate all appropriate redshift effects. Then the mean expected 
flux x a at a point A is ^[dx^/ dri\(dn) , while the expecta- 
tion value of the flux at point A times the flux at point B 
is dn\\\dx B a l ' dn2\,(dn\dn 2 ) ■ When the number of 
small volume elements is large, the double sum over differ- 
ent halos dn\ ^ dn 2 yields an expectation value x^x^ which 
is simply x 2 a , while the sum over dn\ = dn 2 (which we then 
denote simply dn) yields the correlated Poisson term. Thus, 
the correlation function of relative Poisson fluctuations in flux 
is 

«o4EttW' (2l) 

x z a — ' dn dn 

where the points A and B are separated by a distance I. 

To apply this equation, we divide halos into effective four- 
dimensional volumes dVdM, which include three spatial di- 
mensions as well as a distribution in halo mass, and convert 
the sum into a four-dimensional integral. In a homogeneous, 
non-expanding universe, the only difference between dx^/dn 
and dx B a /dn would be the distance factors l/r\ and l/r|. For 
distant sources, r& ~ rs ^> I, we have an integral dV /r 4 at all 
r > I, so the contribution to £p(Z) of all such sources is oc I /I. 
On the other hand, sources near A satisfy r& <C / and ~ I, so 
their contribution is dV /{r\l 2 ), over r& from to ~ /, which 
is again oc l/l. Sources near B contribute equally, and the 
overall correlation function diverges at small separations as 
£p(Z) oc 1 /Z, as expected from the simple scaling argument at 
the end of §|3 

To apply this formalism in the present context, we first write 
the flux as a double integral over the the halo number density 
as a function of mass and volume. Halos in a comoving vol- 
ume dV at r contribute a flux oc 1 jr 2 at the origin, and we 
assume that the fraction of the flux contributed by halos of 
a given mass M is proportional to the gas fraction contained 
in these halos. Thus, we write the contribution dx a from the 
halos in dV as 

, 1 f dn(z') 
dx a = P(z') -= \ M—^-dMdV , (22) 
r Jm dM 

where we have factored out expl icitly the 1/r 2 dependence. 
Using the Sheth & Tormen ( 1999) halo mass function dnjdM 
at z 1 , and dx a calculated as in § the above equation de- 
fines P(z?) as an overall normalization factor that accounts for 
the spectral distribution and redshift effects. Note that the 
mean intensity x a at a point is given by the volume integral of 
equation ( 1221 . which yields the same result as given by equa- 
tion iQJ) . 

An additional subtlety necessitates a small adjustment in 
equation J2 It before we can apply it in the present context. 
This equation, as written, requires that the same halos be seen 
by points A and B at the same redshift z- However, if, e.g., 
point A is closer to some location r, then the point B will see 
the progenitors of the galaxies at r seen by point A. When we 
apply equation d21i . we integrate over the halo distribution 
seen by the closer point A, and note that the total emission 
rate of its progenitors at the higher redshift z' B = z'(?b) is on 
average some fraction of the emission rate of the later parent 
galaxy at z' A = z'(ot)- An actual correlation function measured 
over a significantly large volume will indeed average over the 
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progenitor distribution. This average fraction may depend on 
parent mass and on the two redshifts. First, we neglect the de- 
pendence on parent mass, i.e., we consider a mass-averaged 
fraction. This is likely to yield an accurate result at high red- 
shift, since when galaxies are rare they are all found in a nar- 
row range of halo masses just above the threshold mass for 
star formation. We then fix the redshift dependence of this 
factor so that the total gas fraction in progenitors at z' B is con- 
sistent with the total fraction calculated in the normal way 
from the halo mass function at z' A . Thus the dimensionless 
correlation function is 



£p(0 = 4- 



dV 



dn(z' A ) dMM2 P(z' A )P(z' B )F s M' B ) 



M 



dM 



^gal(^) 



(23) 

where dV =d s rA, and we integrate over half the volume where 
fA < r B, with the factor of 2 accounting for the contribution of 
sources that are closer to B. 

We have assumed in calculating r B that the progenitors of 
a galaxy are at the same location as the parent. The actual 
displacements are at most of order the Lagrangian size of a 
typical galactic halo, ^0.1 comoving Mpc, which is negligi- 
ble when r B 3> 0.1 Mpc, while if we consider small values of 
r B then the displacements resulting from gravitational infall 
velocities (< c/1000) during the time interval \r B — rj^/c are 
still a small fraction of the distance r B -rA- 

The correlation function £p(7) drops sharply at large dis- 
tances, and in particular, when I is more than twice the dis- 
tance corresponding to z max (2), no halo can contribute flux to 
both points A and B and thus £p(l) = 0. At small distances, 
£p(/) oc 1 // as in a non-expanding universe, but even for I — > 
the proportionality factor is affe cted by dist ant sources since 
they dominate the flux [compare Zuo ( 1992)]. Figure ^illus- 
trates this correlation function in several cases. The maximum 
Lya horizon in this case, corresponding to the redshift z ma x(2), 
is 284 comoving Mpc. When we include only the stellar ra- 
diation emitted up to Ly/3, £p(l) begins to drop below the 1// 
asymptotic form at separations / > 200 Mpc and sharply drops 
to zero at twice the Lya horizon. However, when more ener- 
getic radiation is included, £p(7) drops below the 1 / 1 power 
law at much smaller I, corresponding to the smaller effective 
source horizons of the energetic photons, thus making the cut- 
off at I ~ 200 Mpc much rounder and less prominent. 

In this Figure, the star formation parameters are fixed to 
produce the Lya coupling transition (i.e., x tot = 1) at z = 20. 
In our standard case of full Lyman-band emission and Pop 
III stars, with only atomic cooling occurring in galactic halos, 
the required star formation efficiency is relatively low, 0.22%, 
and / esc = 0.3 yields an ionized fraction of 1% at z = 20, with 
complete reionization at z = 9. Although here we have calcu- 
lated £p(/), in the rest of the paper we Fourier transform the 
Poisson-induced correlation function and find the correspond- 
ing power spectrum. 

For the case of a mass-dependent star formation efficiency, 
we insert the appropriate weighting factor into the various re- 
lations given in this section. Note that if star formation in in- 
dividual galaxies actually proceeded in short, random bursts 
with a duty cycle f|« 1, then the higher rate of star forma- 
tion during bursts (by a factor l/r/) would counterbalance the 
fraction of time of burst activity (77), and the overall correla- 
tion function would be unchanged for both of the sources that 
we consider for fluctuations in flux. 

We insert two additional cutoffs that are expected on 
small scales. The first is the cutoff due to baryonic 
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FIG. 3 . — The Poisson-induced dimensionless correlation function of 21cm 
brightness fluctuations versus comoving separation. We assume that galaxies 
formed via atomic cooling in halos at z = 20 with a fixed star formation effi- 
ciency in all halos, set to produce the coupling transition at this redshift (so 
that x t oi = 1 )■ We also assume that the IGM gas has been cooling adiabatically 
down to this redshift. In each set of solid curves we include, from bottom to 
top at / = 1 Mpc, stellar radiation emitted up to Ly/3, Ly<5, or full Lyman-band 
emission (see §|5J, assuming Pop III stars. Also shown for comparison is the 
case of full Lyman-band emission with Pop II stars (dotted line). We also 
show the separation equal to twice the maximum Lya horizon (dashed line). 



pressure, which affects the density fluctuations and has 
an approxima t e sha pe (1 + k 2 R\)~ 2 ( Sh apiro et alJ 1 1994 
Gnedi n & Huil 119981) where the filtering scale is Rp = 

2.9 (M F /10 6 M Q ) 1/3 (O m /! 2 /0.14)" 1/3 kpc, in terms of the fil- 
tering mass Mp which is a weighted time average of the Jeans 
mass(|Gnedin & Hui 1998). We calculate M? appropriately for 
the two cases we consider, an adiabatically cooling IGM or 
a pre-heated IGM. The second cutoff that we include results 
from the thermal width of the 21cm line which also smooths 
out the fluctuations on small scales. We approximate this ef- 
fect as a Gaussian cutoff, e~ k * T , where the spatial scale corre- 
sponding to the one-dimensional thermal velocity distribution 
is 



R T = 7.0 



( m P \ 
\l.22) 



-1/2 



100K 



1/2 



1+z 
10 



-1/2 



0.14 



-1/2 



kpc. 



(24) 

and m p is the mean mass per particle of neutral primordial gas 
in units of the proton mass. 

For the same parameters as in Figure|5J Figure |4] shows the 
/i 2 angular component of the observable power spectra, given 

by 



PAk) = 2P s (k) 



Xtot 



(25) 



and P lln _,5(A:) which equals the Poisson power spectrum except 
for a factor of (x Q /x to t) 2 . In addition to these quantities, P^o 
yields the baryonic density power spectrum at each redshift, 
and its normalization (measured in mK) thus yields the mean 
brightness temperature 7},. Note that the actual plotted quan- 
tity is \Ti ! \[k i P(k)/(2Tr 2 )] 1 ^ 2 , i.e., the square root of the contri- 
bution per log k of power on the scale k to the variance of the 
brightness temperature. The figure shows that (3 can be deter- 
mined from Ppi/P^i on small scales (k > 1 Mpc" 1 ), making 
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FIG. 4. — Power spectra of 21cm brightness fluctuations versus comoving 
wavenumber. Shown are P 2 (top panel), which contains the contribution 
of the density-induced fluctuations in flux, and P un s (bottom panel), which 
is solely due to the Poisson-induced fluctuations. We assume that galaxies 
formed via atomic cooling in halos at z = 20 with a fixed star formation effi- 
ciency in all halos, set to produce the coupling transition at this redshift (so 
that Jttot = 1). We also assume that the IGM gas has been cooling adiabatically 
down to this redshift. In each set of solid curves we include, from bottom to 
top at k = 0. 1 Mpc~' , stellar radiation emitted up to Ly/3, Ly<5, or full Lyman- 
band emission (see §|5}, assuming Pop III stars. Also shown for comparison 
is the case of full Lyman-band emission with Pop II stars (dotted line). We 
also show 2(3Pg (dashed curve). 



the two cases (adiabatic: j3 < 1/3 or pre-heated: j3 f» 1) easily 
distinguishable. As a result, the mean brightness temperature 
determines x tol in these two limits [see eq. ©]. Measurements 
at k > 100 Mpc -1 can also independently probe 7i because 
of the smoothing effects of the gas pressure and the thermal 
width of the 21cm line. 

Note that the value of x a (which dominates jc tot in the cases 
shown in Figure |3 is interesting in itself, since it fixes the 
product of the fraction of gas in galaxies, times the typical 
efficiency with which gas is turned to stars and the stars pro- 
duce ultra-violet photons inside these galaxies. The Figure 
also shows that the power spectra are strongly sensitive to 
the radial distribution of flux seen by the gas, which in turns 
probes the contribution of photons from the higher-level Ly- 
man series, a contribution which depends on the stellar IMF 
and the spectrum of the sources. In particular, since £/>(0 
drops to zero at large /, the corresponding Poisson-induced 
power spectrum approaches a constant at small k, and the plot- 
ted quantity is oc 

X-ray heating of the cosmic gas above the CMB tempera- 
ture is likely to occur well before the universe is ionized be- 
cause it requires an amount of energy of only ~5x 10~ 3 [(1 + 
z)/20] eV per baryon, while reionization requires ~ 13.6/"' 
eV per hydrogen atom (or even more if recombinations are 
included). The X-rays may be produced by stellar winds, 
supernovae, X-ray binaries, quasa rs, or collapsed halos fse e 
Madau et al] 11997ft and §2.3.2 in IWvithe & Loebl J2004a)l. 
As long as more than ~ 0. 1 % of the energy radiated by galax- 
ies above the Lyman-limit couples to the IGM through X-ray 
scattering, 7^ exceeds T 7 before the completion of reioniza- 
tion. It is possible, however, that Lya coupling between T s 
and Tk is achieved before X-rays heat Tk above T 1 . 



We note that a significant gas fraction may already lie 
in virialized minihalos at high redshift. Although this gas 
may be dispersed due to internal feedb ack from star forma- 
tion due to molecu lar hydrogen cooling \OYi & Haimanl2003T: 
iRicotti et al.ll2002l) . if stars cannot form then th e dense ga s 
minihalos produce 21cm fluctuations (Iliev et al. 2002, 2003). 
Since the gas content is strong ly suppressed in halos below the 
filtering mass JGnedinll2 000). we estimate the maximum gas 
fraction in minihalos as lying between the filtering mass Mp 
and the minimum mass needed for atomic cooling. At z = 20, 
only around 0.6% of the gas is in minihalos in the mass range 
8 x 10 5 - 3 x 10 7 M Q , compared to 11% in minihalos in the 
range 4 x 10 s - 1 x 10 8 M Q at z = 8. Thus, in the redshift 
range that we focus on in this paper, only a minor fraction of 
the gas is in virialized regions and so we consider the density 
and temperature distributions only in uncollapsed regions. 

7. STUDYING THE EARLIEST GALAXIES 

Observing the enhanced 21cm fluctuations would con- 
stitute a direct detection of the collective emission of the 
earliest generation of galaxies, long before reionization. 
Such a detection is potentially feasible within a few years 
(http://space.mit.edu/eor-workshop/), while the number of 
stars in such an early galaxy may be smaller than in a 
present-day globular cluster, making these galaxies at z ~ 20- 
30 undetectable individually even with the proposed James 
Webb Space Telescope (http://www.jwst.nasa.gov/) planned 
for launch in the next decade. 

In order to illustrate how the properties of the first galax- 
ies may be probed with 21cm fluctuations, we predict the ob- 
served power spectra at a sequence of redshifts, for various 
choices of the parameters of galactic halos and of IGM heat- 
ing. Figure |3 shows a reasonable lower limit on the number 
of halos in the universe. Here we maximize the dominance 
of massive halos by assuming that: (i) atomic cooling is re- 
quired for efficient star formation, and ( ii) the star formation 
efficiency de clines with circular velocity according to local 
observations JDekel & Wool2002t lkauffman et al. 2003J, em- 
phasizing the role of high-mass halos. The star formation ef- 
ficiency is set to 12% above V c = 180 km s" 1 and oc V 2 below 
this value. In this case, reionization completes at z = 11.7. In 
Figure [5] the threshold parameter x lot (going from high red- 
shift to low) is 0.070, 0.037, 0.12, 1.0, and 6.2. The corre- 
sponding value of /3 is 2.1, 2.0, 0.36, 0.22, and 0.25, respec- 
tively. We find at z = 20 a filtering mass M F = 8.0 x 10 5 M Q , 
Rf = 2.7 kpc, and Rj = 1.4 kpc. Collisions dominate the 
Lya coupling at z > 30 (contributing x c = 0.070, 0.029, and 
0.008 at the three highest redshifts shown), producing a high 
/3 and suppressing the effect of the fluctuations in Lya flux 
which are largest at the highest redshifts. The absolute signal 
in mK tends to rise with time due to the growth of density 
fluctuations (which affects P^i) and the declining gas tem- 
perature which increases the overall 21cm absorption in the 
adiabatic case. However, the increasing number density of 
sources with time eventually suppresses the Poisson fluctua- 
tions. Also, the saturation of the Lya coupling (i.e., x tot 3> 1) 
reduces the sensitivity of the 21cm fluctuations to fluctuations 
in the Lya flux; in the figure, by z = 15 the power spectrum 
P^iik) has essentially reverted back to the shape of Ps{k). 
The figure does not include the signature of ionized bubbles, 
which becomes significant only at later times when the ion- 
ized fraction approaches unity. 

Figure [6] shows a reasonable upper limit on the number of 
halos in the universe. In this case we maximize the contri- 
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FIG. 5. — Power spectra of 21cm brightness fluctuations versus comoving 
wavenumber. Shown are P 2 (top panel), which contains the contribution 
of the density-induced fluctuations in flux, and P un s (bottom panel), which 
is solely due to the Poisson-induced fluctuations. We assume atomic cool- 
ing with an efficiency that varies according to local observations (see text), 
normalized so that „t t oi = 1 at z = 20. Redshifts shown are z = 35 (dotted), 
30 (short-dashed), 25 (long-dashed), 20 (upper solid), and 15 (dot-dashed). 
We also show 2f3Pg at z = 20 (top panel, lower solid curve). We assume full 
Lyman-band emission (see §|5} and Pop III stars. 
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FIG. 6. — Power spectra of 21cm brightness fluctuations versus comoving 
wavenumber. Shown are P^i (top panel), which contains the contribution of 
the density-induced fluctuations in flux, and P un s (bottom panel), which is 
solely due to the Poisson-induced fluctuations. We assume H2 cooling with 
a fixed efficiency set so that x to t = 1 at z = 20. Redshifts shown are z = 35 
(dotted), 30 (short-dashed), 25 (long-dashed), 20 (upper solid), and 15 (dot- 
dashed). Note the different axes ranges compared to Figure|5] We also show 
2f3P$ at z = 20 (top panel, lower solid curve). We assume full Lyman-band 
emission (see §|5) and Pop III stars. 



bution of numerous low-mass halos by assuming that only 
molecular cooling is required for efficient star formation, with 
a fixed star formation efficiency (/» = 0.03%). With these 
parameters, the universe would not complete reionization by 
redshift 6 without a significant rise in /* at z > 6. Figure |6] 
shows some of the same trends that are apparent in the pre- 
vious figure, although both sources of fluctuations in flux are 
significantly smaller, since there are more halos contributing 
and they are more weakly biased. The threshold parameter x lot 
(going from high redshift to low) is .076, .082, 0.30, 1.0, and 
1.9. As before, collisions dominate at high redshift, contribut- 
ing x c = 0.070, 0.029, and 0.008 at the three highest redshifts 
shown. The values of are 1.9, 0.95, 0.24, 0.22, and 0.25, 
respectively. 

In FigureQwe consider the same halo population as in Fig- 
ure |5| but with the 21cm fluctuations modified due to an as- 
sumed pre-heating of the IGM. We set the star formation ef- 
ficiency to 6.9% above V c = 180 km s" 1 and oc V c 2 below this 
threshold, so that x tot = 1 at z = 20. The values of x a at the 
various redshifts are 0.0003, 0.006, 0.08, 0.79, and 5.3, while 
x c = 0.61, 0.46, 0.32, 0.21, and 0.12, respectively. In this case, 
collisions are still significant at the lower redshifts because of 
the increased collision rate at high temperature. The value of 
[3 is 1.6, 1.7, 1.6, 1.1, and 1.0, respectively, always > 1 since 
the gas is uniformly heated and no longer changes its temper- 
ature adiabatically. In this case M F = 2.4 x 10^M Q , R F = 3.9 
kpc, and Rj = 10.8 kpc at z = 20. Comparing Figures [5] and 
we see that the strong differences between these two cases 
in both the mean 7j and in the value of (3 (which affects P^) 
allow for a consistency check on the thermal state of the IGM. 
Note that in the cases shown in Figures [5] and the Poisson 
fluctuations diverge on small scales, in which case our lin- 
ear analysis based on equation @ represents an estimate that 
is subject to strong non-linear corrections. Also, on scales 
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FIG. 7. — Power spectra of 21cm brightness fluctuations versus comoving 
wavenumber. Same as Figurel5l except that we assume pre-heating to 500 K 
by X-ray sources. Note the different axes ranges compared to Figurel5l 



< 0.1 comoving Mpc, which corresponds to the initial region 
that collapses to become a typical halo at z ~ 20, our calcu- 
lation is subject to corrections due to non-linear gravitational 
collapse. 

From the ratio P^i/P^t at large scales, it is possible to mea- 
sure W(0) [see eq. J25H . which yields directly the average bias 
factor of the emitting sources, given by equation J20l > with the 
curly brackets becoming simply {1 + b(z')} in this limit. The 
relative bias of galaxies is on large scales predicted to depend 
simply on the rarity (# of sigma) of the underlying Gaussian 
density fluctuations that led to their formation. Thus, it is sen- 
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sitive to the typical mass of the dark matter halos that host 
galaxies at these redshifts. This large-scale enhancement fac- 
tor can be easily seen by comparing the upper solid line to the 
dashed line in Figure [4] and the pair of z = 20 lines in each 
of Figures|5j|6j and0 The enhancement factor shown in the 
figure (which equals the square root of the ratio between P„i 
and 2f3Ps) is in each case 4.8, 5.1, 4.0, and 2.2, respectively. 
The corresponding value of W(Q) is 9.8, 11.2, 6.5, and 11.2, 
respectively. 

8. SUMMARY 

We have shown that the 21cm flux fluctuations from in- 
tergalactic hydrogen are strongly enhanced due to the radi- 
ation emitted by the first galaxies at redshifts ~ 20-30. The 
earliest galaxies emitted photons with energies between the 
Lya and Lyman-limit transitions of hydrogen, to which the 
neutral universe was transparent except at the Lyman series 
resonances. As these photons redshifted into a Lyman-series 
resonance and then degraded into a Lya photon, they cou- 
pled the spin temperature T s of the 21cm transition of hydro- 
gen to the gas kinetic temperature 7i, allowing it to deviate 
from the microwave background temperature, T 1 . Figures 4— 
7 show that the fluctuations in the radiation emitted by the 
first galaxies enhanced the overall power-spectrum of 21cm 
flux fluctuations by at least a factor of a few during the pe- 
riod when Lya coupling was established. The enhancement 
is sourced by inhomogeneities in the biased distribution of 
galaxies (as those represent rare peaks in the density field) 
as well as by Poisson noise in the number of galaxies. The 
effects of galaxies are more apparent when the power spec- 
trum is split into angular components relative to the line of 
sight (with the anisotropy sourced by peculiar velocities); in 
particular, the component P un _,5 shown in the lower panels of 
Figures 4-7 (eq. II It would vanish in the absence of Poisson 
fluctuations in the distribution of galaxies. For the density- 
induced fluctuations (upper panels), the enhancement factor is 
roughly constant on scales > 10 comoving Mpc. Because of 
the rapid evolution in the abundance of galaxies at early cos- 
mic times, and the atomic physics of resonant scattering by 
hydrogen, the fluctuations in the Lya intensity are dominated 
by the contribution from sources within distances ~ 1-100 
comoving Mpc, well within the source horizon out to which a 
Ly/3 photon could redshift into a Lya photon (Fig.|2ji. 

Several experiments are currently being con- 
structed or designed to detect the 21cm fluctua- 
tions prior to reionization. For an overview, see 
http://space.mit.edu/eor-workshop/, for individual ex- 
periments, see LOFAR (http://www.lofar.org/), MWA 
(http://web.haystack.mit.edu/arrays/MWA/index.html), 



PAST (http://xxx.lanl.gov/abs/astro-ph/0404083), and the 
future SKA (http://www.skatelescope.org/). Although low- 
frequency foregrounds are much brighter than the 21cm 
signal, they are not expected to include sharp spectral fea- 
tures, unlike the 21cm maps that should decorrelate over small 
shifts in frequency (corresponding to slicing the universe 
at different redshifts). Consequently, the prosp ects for ex 
tractio n of the 21cm sign a l are p r omising ((Gnedin & Shave i 
2004; Zaldarr iaga etall l2004t iMorales & Hewittl 1200- 
Moralea l2004t ISantos et alJl2004l) . The signal-to-noise ratio 
in these experiments is controlled by the sky brightness and 
grows as the square-root of the integration time. Hence, the 
increase in the amplitude of 21cm fluctuations that we predict 
has important practical implications for observing the 21cm 
signal from the corresponding redshifts, since the required 
integration time is reduced by the square of the enhancement 
factor. 

The enhanced level of fluctuations in the 21cm flux lasted 
only until the Lya coupling of T s and 7* was saturated. Both 
before and after the Lya coupling epoch, the enhancement 
(compared to the case of a uniform Lya flux) had a lower am- 
plitude (see, e.g., Figure [5] where Lya coupling was achieved 
at z ~ 20). Detection of the enhancement factor over the 
available redshift interval (i.e., the corresponding observed 
frequency band) can be used to infer the abundance and 
characteristic mass of the earliest galaxies (§7) without any 
such galaxy being observed individually. The existence of 
Lya coupling can also be used to determine whether 7^ is 
larger than T 7 in this redshift interval. If the 21cm signal ap- 
pears in emission rather than in absorpti on fwhere abso rption 
is pre sent in the interval 30 < z < 200 (lLoeb & Zaldarriaga 
2004)] then 7* > T~, most likely due to heating from X-rays 
emitted by the same galaxies. Altogether, the 21cm fluctua- 
tions can be used to probe the level of UV and X-ray radiation 
emitted by the earliest galaxies, long before reionization was 
complete. The milestones of Lya coupling between T s and 
Tk and X-ray heating of 7* above T 7 provide important pre- 
cursors for the radiative feedback of galaxies on the IGM that 
culminated with the completion of reionization. 
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